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FOR A SPIN 1 PARTICLE IN EXTERNAL COULOMB FIELD 
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Quantum-mechanical system - spin 1 particle in external Coulomb field is studied on 



{^JQ. the base of the matrix Duffin - Kemmer - Petiau formalism with the use of the tetrad 

technique. Separation of the variables is performed with the help of Wigner functions 



D-mai^i 0, 0); a ~ —1, 0, +1; j and m stand for quantum numbers determining the square 
and the third projection of the total angular momentum of the vector particle. With the help 
of parity operator, the radial 10-equation system is divided into two subsystem of 4 and 6 
Qh' equations that correspond to parity P = (—1)^+^ and P — (—1)-' respectively. The system of 



4 equation is reduced to a second order differential equation which coincides with that arising 
in the case of a scalar particle in Coulomb potential. It is shown that the 6-equation system 
reduces to two different second order differential equations for a "main" function. One 
main equation reduces to a confluent Heun equation and provides us with energy spectrum. 
QQ ' Another main equation is a more complex one, and any solutions for it arc not constructed 

oo 
vn 

od 

O ; I. INTRODUCTION AND SEPARATION OF THE VARIABLES 

^ Many years ago, a very peculiar behavior of a spin 1 particle in presence of the external Coulomb 

^ ' field was noticed by I.E. Tamm As far as we known the whole situation with this system stays 



much the same. In the present paper, we start examining the problem on the base of the matrix 

5 [with the use of the tetrad generally covariant 



Duffin - Kemmer - Petiau formalism [^, Q], [J, 
technique (for more details and references see in 



16-13, 



it turns out to be more convenient 



than a common Proca tensor approach Q| ) • Choosing a diagonal spherical tetrad according to 

dS^ = dt^ - dr^ - r'^ide^ + sin^ Od^'^) , 
efo) = (1,0,0,0) , ef3) = (0,1,0,0) , 
e?i) = (0,0,-,0) , 6%) = (1,0,0,— ^) , (1) 
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we reduce the main D-K-P stationary equation to the form 



^{x) = , 



(2) 



where e = E/{ch),a = e'^/{ch),M = mc/Ti and Se^^ stands for an angular operator (its form 
means that we have here a generahzed Schrodinger - PauU basis 

2 id + i j^^ cos 



sm ( 



Krohkowski and Turski 
n 



Spherical waves with (j, m) quantum numbers should be constructed within the following general 
substitution (we adhere notation developed in Red'kov [16-22, 301; before similar techniques was 

M], Turski [15] ; many years ago such a tetrad 



applied by Dray jl^ 

basis was used by Schrodinger [8'] and Pauli ^\ when looking at the problem of single- valuedness 
of wave functions in quantum theory ~ then the case of spin 5 = 1/2 particle was specified) 

^{x) = { ^q{x), $(x), E{x),H{x) } , 



M^) = e-^^*/^ Mr) Do , Hx) 



-iEt/n 



E{x) = e-^^^'^ 



Ei{r) D^i 
E2{r) Do 
Ez{r) D+i 



H{x) = e"^-^*/^ 



short notation for Wigner functions ^0[| is used: D(^ = D\ 



recurrent formulas 



accordance with Pauli approach the quantum number j takes values 0,1,2, ... With the help of 



$i(r) 


D-i 


^2{r) 


Do 


Mr) 


D+i 


Hi{r] 


D-i 


H2{r] 


Do 




D+i 




0) , 



(3) 



ul (where u = + 1) , a = V(j-l)(j + 2)) 



]^ — COS 1 

de D^i = - {aD^2-i^ Do) , D^i = - {a D^2 + ^ Do) , 

I sm d Z 



de Do = ]^{v D_i - V D+i 



"fn ^ 1 . „ „ 
Do = -{vD_^^vD^^ 
sm ij 1 



jj^ _|_ (3Qg Q \ 

9e D+i = - {u Do - a D+2 ) , ^ — D+i = - { u Do + a D+2 ) , 

2 sm & 2 



(4) 



after simple algebraic calculation we arrive at the radial equations (for clarity corresponding Proca 
tensor equations are written down as well) 
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{-f + - E2 

ar r 



{El + E3 



a, , d 1, u 

+i{e + -) El + Hi+ i- H2 

r ar r r 

+i{e + ^) E2 - i"- {Hi - H3) = M$2 , 

a, , d 1^ V 

+i(e ^-)E3 - i {— + -) H3 - i- H2 = M <^>3 ; 

r dr r r 



ab 



i{e + -Hi + - ^0 - MEi 

r r 



a. 



d 



dr 



$0 - ME2 = , 



-i(e + -) ^2 
r 

-i{e +^)'^3 + ^ «>o - ME3 = , 

-i (-^ + -)$!- i- <^2 - MHi = , 
dr r r 

+i- ($1 - $3) - MH2 = , 
r 

+i {4- + -)^3 + - MH3 = , 
dr r r 



(5) 



(6) 



where v = \/ j{j + l)/2 (note the factor 1 / y/2) . 

Concurrently with J^, J3 let us diagonalize an operator of spatial inversion 11. After transition 
to spherical tetrad basis, and also to cyclic representation for D-K-P matrices /3", for this discrete 
operator we get 



n 



100 

Hs 



-U3 



P 



n. 



0-1 
0-1 
-100 



(7) 



Eigen-value equation n\I' = P ^ results in two different in parity states 



P=(-l) 



^>o = , $3 = -$i , $2 = 0, 

E3 = —El , £'2 = 0, H3 = Hi ; 



(8) 



P = {-ly , ^.3 = $1 , = +Ei , H3 = -Hi, H2 = 



(9) 



Correspondingly, eqs. ([5]) - ([6]) give 4 and 6 equations 
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+ i(e + -) Ei + i{4- + + ^-^^2 = M$i , 

r ar r r 

-i(e + -)^i = MEi, + = MHi , 2i-^i = MH2 , 

r dr r r 

excluding Ei,Hi,H2 we get a second order differential equation for 



<f 2d , Q,9 , ,9 if? + 1) 
ar^ r dr r 



«>i = 



(10) 



(11) 



which coincides with that arising in the case of a scalar particle in Coulomb potential (it is the 
known fact that was noted in Tamm's first paper jlj]. Its solution is well known and provides us 
with the following energy spectrum (in usual units) 



E 



nic 



1 



For states with parity P = (—I)-' we have the system 



N = n+- + .y{j + 1/2)2 _ ^2 _ 



(12) 



-ly, 



(4- + -)E2 + 2-Ei + Af$o = , 
dr r r 

+i{e + -) Ei+ i{— + -)Hi - M$i = , 
r dr r 

+i{e + ^)E2 - 2i'^Hi - M$2 = , 

Ct V 

-i{e + -)$! + -$o - MEi = , 

(y. d 
i{e + -)$2 + + ME2 = , 
r dr 

d \. u 
i{-r + -)$i + «-^2 + MHi = . 
dr r r 



(13) 



II. THE CASE OF MINIMAL VALUE j = 



The states with minimal value j = can be treated straightforwardly (also see in [20]). In this 
case, we should start with special substitution of the wave function 

<i>o(x) = e-^^*/'*$o(r) , C{x) = e-^^*/'*C7(r) , Co(x) = e'-^Co (r), 



C{x) = e-'^'l^ 










E2{r) 




<l>(x) = e 



-let 



H{x) = e 





^>2(r) 







(14) 
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The Eg acts on this function as a zero operator, and parity P = (—1)'^+^ = —1. Now the radial 
system is (for ehminating imaginary unit i, we use slightly different variables: $o = fO: —i^i = 



H2 = 0, -(4- + -)E2 = Mipo , 
dr r 

(y ad 

(e + —)E2 = Mip2 , (e H — )ip2 z-ipo = ME2 . 

r r dr 



(15) 



Prom whence it follows a second order equation for E2 



flj.2 J. J.2 



E2 = 



(16) 



With substitution E2{r) = r ^f{r), we get 

dr"^ 



(17) 



In dimensionless variables 



rE rn^c* 



X = re = 



E^ 



= A^ 



it reads 



ax^ X x'^ 



With the substitution f{x) = x"e ^^F{x), for F we obtain 



)/ = 



(18) 



xF" + (2a - 2hx)F' + [ ^ "'^ ^ + (6^ + 1 - A^)^ + (2a - 2ah) ]F = . 



X 



Requiring 



1 ± V9 - 4a2 — - , Vm2c4 - ^2 
a = ^ , = ibv — 1 = =b- 



(19) 



2 , - v- E ' 

the choice of upper signs in the formulas provides us with good parameters for bound states, we 
get 



X F" + 2(a - hx) F' + 2(a - a6) F = . 



(20) 



This equation is solved in confluent hypergeometric functions. Let us specify this solution in detail. 
With the use of expansion for F{x) 



F{x) = ^CkX^ , F' = kCk x^-^ , F" = J2 Hk - 1)Q X 



k-2 



k=0 



k=l 



k=2 



6 



we get 

oo oo oo oo 

k{k - l)Ck x^-^ + 2a ^ kCk x^'^ - 26 ^ kCk x^ + 2(a - ab) ^ Ck x^ = , 

k=2 k=l k=l k=0 

from whence it follows 

[ 2aCi + 2(a - ab) ].'r" + [2C2 + 2a 2C2 - 26 Ci + 2(q - ab)Ci ] x + 

00 

+ ^ [ n(n + l)C„+i + 2a(n + - 26raC„ + 2{a - ab)Cn ] a;" = . 

ra=2 

Therefore, we arrive at recurrent formulas 

Ci = -{a-ab) Co = , 
C2 2 (1 + 2a) = 2 [6 - (a - a6)] Ci = 0, n = 2, 3, 4, 
Cn+i (n + 1) (n + 2a) = 2 [n 6 - (a - a6) ] C„ = . (21) 

To get the polynomial we must require 

Cn+1 = ^ [ iV 6 - (a - a6) ] = . 

This quantization rule gives 



a-a6 _ 1 + V9 - 4a2 Vm2c4 - ^2 

— r — = iv , a = , 6 = - 



E 



so that 



2ae - (1 + V9 - 4:a^)Vm'^c'^ - E^ _ ^ _ 
2Vm'^c'^ - E"^ ~ ' 



its solution is 



2 

E = = , 2r = 1 + V9 - 4a2. (22) 

VI + a2/(r + Ar)2 ^ ^ ' 

III. LORENTZ CONDITION IN PRESENCE OF COULOMB POTENTIAL 

As known for a massive spin 1 particle there must exist a generalized Lorentz condition. Let 
us specify it for the problem under consideration. In the Proca tensor equations for the vector 
particle 

^p-Dp^^ = M , = M , (23) 
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where Da = Vq + i (e/c) Aa, let us act by the operator Da on the second equation in ([23]) . we get 

(V. + i -Aa) = ^ Fa^ -^"-^ ■ (24) 

This Lorentz condition should be translated to tetrad form. To this end, instead of and ^'"^ 
one should introduced their tetrad components 

Correspondingly, eq. (j24p takes the form 

(eW" ,1,^ + ^{a). da ) ^(.) + i I A^'^) ^(,) = i ^ FW('') . (25) 

The Coulomb field Aq = e/r , Fro = — e/r^ in the tetrad description looks 

^(0) ^ g(o)o^^ ^ £ ^ ^30 ^ e(3)'^eW° F,o = - ^ . (26) 

Beside, after simple calculation we get 

g(0)a^n COS g e(2)« = n e(3)A^_2 

'° ' '° rsine ' ^ ;° ^ ' ^ ^ • 

The components functions '^(a) ^"^^ ™ ^25p can be related with components of the D-K-P 

column as follows (transition between cyclic and Cartesian representations; c= l/\/2) 



^(0) 




1 































$0 A 


^(1) 







— c 





+c 






















$1 L>_i 









—ic 





—ic 






















$2 A 


^(3) 










1 

























$3 L>+i 


^(0)(1) 
















— c 





+c 













El D^i 


^(0)(2) 
















—ic 





—ic 













E2D0 


^(0)(3) 



















1 
















E3D+1 


^(2){3) 

























— c 





+c 




HiD^i 


^(3)(1) 

























—ic 





—ic 




H2D0 


^(1){2) 




























1 







H3D+1 



We need only $0, *(i), ^(2), ^(3)> ^(o)(3) : 

= e-^^*/^-^ (-<I>i D^i + «>3 D+i) , 

*(2) = e-*""*/^^ (-^1 D-i - $3 D+i) . 
*(0) = e-'^'/^^o Do , ^(3) = e-*^*/^<I>2 , 

*(0)(3) = e-'""'/^ E2 Do . (27) 
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With the help of (|27]), eq. gives 



1 M-lcos9 

V 2 sm 

1 ^ ^ M + lcos6') ^ , 

_ r $3 {de D+i + — D+i) + 

V 2 sm 

,2 d , a 

+Do (-- - ie c^o - ^ ^2) = i i?2^o 



Now, by taking into account the recurrent formulas 

M - 1 cos I 



111 



de D^i — D^i = -^/(J + Tjj Dq , 

M + 1 cos 



de D+i D+i = -VU+Tj] Do , 

sm 

we arrive at a the Lorentz condition in radial form 

--^0- (|:+')^2-^('l>i + <&3) = ^i^2. (28) 

For states with parity P = (—1)-''''^, this condition is satisfied identically. For states with parity 
P = {—ly it gives 

/ 2, 2i/ ^ ia ^ ,^„, 

_,,.f„_ + (29) 



From eq. (I29p a very important relationship can be established. To this end , from eq. (I29p let 
us exclude $2 with the help of the third equation in (fT3]) 



, , , , a. , d 2 , ^ 2iz^ , d 1 , ^ 2z^M , a ^ 

I eM «>o + z(e + -){-r + -)^2 (:r + + $1 = ^ ^2 • 

r ar r r dr r r Iv^ 

Let us transform the second and third terms with the help of 1-st and 2-nd equations in (fT3|) - it 
results in 

E2 = -2Afr $0 • (30) 



IV. THE MAIN FUNCTION OF THE FIRST TYPE, REDUCING THE PROBLEM TO 

THE CONFLUENT HEUN EQUATION 



Let us examine eqs. (|T3]l . now with the use od the Lorentz condition (129)1 and its consequence 
30]) . From the 1-st eq. in p^ . we produce 

A4r d 
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With the help of (j30]) . the fourth eq. in p3|) gives 
$1 



e + ajr \r 2v 
With the help of ([30]), the fifth eq. in ([13]) gives 

$9 



2MV $0, 



$0 • (32) 



^ / , , ■ , ^u, (33) 

The sixth eq. in p3|) provides us with the representation for Hi in term of by nieans of 
second order differential operator 

r 



MH 



, d 1, 
(31 + 



2MV 



$0. (34) 



'dr re + o/ryr 2z^ u dr J re + a/r\dr 

Take note that from the third eq. in p3p one can obtain another representation for Hi that uses 
only the first order operator 

Mr 1 



Hi 



d_ 

dr 



+ 2r [(e+-)2-M2] 



2v e + a/r 

These two representations for Hi must be consistent with each other. 
Now, turning to the Lorentz condition 



.d 2, 2v ^ . a 

one can readily derive a second order differential equation for ^q: 



(35) 



dr^ r \ e + a/r J dr 



a 



+ - ^ - 3M^ + 2M' 



2v' 



e + a/r 



^n = 0. 



(36) 



The function $o will be termed as main function, 5 remaining ones $i, ^>2) -£'1; -f^i are expressed 
through it 



Eo = -2Mr $ 







Mr , d 

Ei = —{f> + 2r—)^o, 
2v dr 



$1 



e + ajr \r 2v 



V 5M2 r^M^ d , 

r 1 $0 



V dr 



$2 



e + Q /r \ dr 



^ - 2MV ) $0 



1 



Mr 



2v e + a/r 



± + 2r[{e + -f-M^ 
dr r 
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(37) 



Changing the variable 



r < 



a 



a 



X ; 



eq. (j36|) is reduced to the form 

1 



dx"^ ' \x X — 1 J dx 
where dimensionless parameters were used 



A2 = a2 a2 



" + ( - , + 1 $0 = 



X - 1 



mc 

A = ^>1. 



Let us consider behavior of the main function near the point x = 0: 

d^^o ^ 3 d^o + 2i^^ ^ 
dx^ X dx x"^ ' 

$0 ~ const x^, A{A - 1) + 3^ - - 2zy2 = , 



-1- \/l + a2 + 2z^2 ^ A = -1 + \/1 + q2 + 2i/2 



(38) 



(39) 



to bound states there correspond positive values of A. In the region near x = +c«, the main 
equation gives 



dx^ X dx 



+ - -T^ + (a^ - A^) $0 = , 



to bound states there correspond solutions vanishing at infinity. 

Now, let us introduce substitution ^o{x) = x^e^^ f{x), eq. (j38|) gives 



dx"^ 



+ 



1 — X 



— + 

dx 



. 2 a2 , 2AB + A + 3B 



+ 



A{A - 1) + 3A - q2 - 2i/2 A + B- 2A2 



x^ 



With restrictions on j4 and B: 



1 - X 



/(x) = 



(40) 



(41) 



A{A - 1) + 3A - q2 - 2z^2 = 



-1+ ^l + 2z^2 + a2. 



5^ + - A^ = 



B = +VA2 - a2 , 



(42) 



eq. (j41|) takes the form 



dx2 



+ 



2B + 



2^ + 3 1 



X — 1 



+ 



2AS + ^ + 35 A + S-2A21 
X 1 — X 



— + 

dx 

f = o. 



(43) 



It can be recognized as the confluent Heun's equation 

J = j{a,o,c,d,h; z) , + { ^ ^ 



dx"^ \ X X — 1 J dx 
[-2d + a{-b - c - 2)]x + a(l + b) + b{-l - c) - c - 2h 



2x{x - 1) 

with parameters given by 



f = 



+2VA23a2, b = +2V1 + 2z/2 + a2 ^ 
c = -2 , d = 2A2 , h = +2 . 



The known condition for polynomial solutions is 

, b+c+2 
d = —a I n + 

it gives the following quantization rule 



{n + Vl + 2i/2 + 0,2)2 



A2 -a2 
Its physical solution is 

2a 



iV2 _ ^N'^ _ 4a2jY2 ■ 
When increases to infinity, we get 

A —7- 00, i/ = m c , ^ m c — ^ = m 

Af2 - Ar2Vl - 4a2Ar-2 2a2 

To obtain a non-relativistic limit, one must impose special restriction 

N = n+ \/l + 2z^2 _^ Q,2 = \/l + j(i + 1) + a2 ^ n + j 

which correlates with the known non-relativistic procedure 

a 

m + eH — wm + e. 
r 

Taking into account eq. ()50p . one can we derive 



1 4a2 , r, . ,^ 



that is 



£;2 = ^2^4 _ (1 + y 1 _ _) ^ ^2^4(^ _ ^) A^ « n + j , 



2 



Thus, the non-relativistic energy levels are given by 

22 4 
, a mc me 

E = 



2N^ h'^m ' 

which coincides with the known exact result (see below). 
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V. THE MAIN RADIAL FUNCTION OF THE SECOND TYPE, REDUCING THE 
PROBLEM TO ANOTHER SECOND ORDER DIFFERENTIAL EQUATION 



Prom general considerations, we may expect two linearly independent solutions for 6-equation 
system for state with parity P = (— I)-'. The above equation (|36p provides us with only one class 
of these, what make us look for another class (may with with some different main function). 

In this connection, let us turn back to eq. ()35p multiplied by — M and compare it with eq. (|34p . 
from that it follows a second order differential equation <I>o, different from (I36p : 

72, 



+ - 



+ 



a 



r(e + a/r) dr 



+ 

la(-15 + 4z^2 
2 r2(er + a) 



,.2 2e^a 



2q2 



er + a r^M'^{er + a) 
e(-5 + 2z^2 _ 3q,2 



r(er + a) 



0. 



In the variable x, it reads 
— -16 



X 



d 



+ 



I J dx 

(-15 + 4z/2 



(1 - X^)a' 



2a' 



+ 



2a' 



X — 1 ' a'^X''x^{x 
(5 + 2Z/2 -Sa^)"! 



2x2(x - 1) 

By means of the coordinate transformation y - 



x{x — 1) 
x~^, eq. becomes 



. 



d^ 



+ 



4y-3 d^>c 



y(i - y) dy 



+ 



:i5 



2a' 



+ 



v'^y 



(1 - X^)a^ 

2/4 2/3(1 -y) ' a2A2(l-2/) 

- 4zy2 + 2a2) (5 + 2zy2 _ 3^2 " 







. 



(53) 



(54) 



(55) 



22/(1-2/) 2/2(1-2/) 
Both differential equations, (I54p and ()55p . are very complex. We might expect that they can 
describe some third class of solutions, however any proof of this does not exist now. 



VI. NON-RELATIVISTIC LIMIT, EXACT ENERGY SPECTRUM 



to 
in 



In^treating this point we will use results on non-relativistic limit for a vector particle according 
22 L general treatment of the problem of the wave equations for arbitrary spin particle see 



21 



24 



25(; also see in [26 1 



First, let us consider the simpler system (I56p for states with parity P = (— l)-'"''"^: 

p = {-iy+\ 



+i{e + -)Ei + + -)Hi + i-H2 = M$i , 
r dr r r 



a. 



-i(e + -)$i = MEi 
r 



+ -)^i = MHi , 2i-<^i = MH2 . 
dr r r 



(56) 
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Here the Hi,H2 represent non-dynamical variables, excluding them we obtain 



I d 1, 

M ^dr r- 



$1 = M$i 



a 



i(e + -)$i = MEi 



(57) 



Now we should make special substitution, introducing a big and small constituents ( Bi(r) and 
Mi(r) respectively) 



$1 = Si + Ml , iEi = Bi- Ml 



(58) 



correspondingly eqs. (j57p take the form 

(e + ^)(^B, - Ml) + ^(|- + \)\B, + Ml) - 1^(51 + Ml) = M(i3i + Ml) 

(e+«)(i^i + Mi) = M(i^i-Mi) 
r 



Summing and subtracting these two we get 

^ ^ ''^ M ^ dr^ ^ """^ Mr^ 

a^,, l,d 1^9, „ 2z.2 



(Si + Mi) = 2MSi , 
(Si + Mi) = 2MMi . 



Now we should separate a rest energy by a formal change e =^ e + M; which results in 



r M ar r 



(Si + Mi) = 0, 



-2{e + ^)Mi + ^ (|: + ^)^ (i?i + ^i) - ^ (^1 + M,) = AMM, . 
Thus, we produce equation for a big Bi{r) and small Mi(r) components 



(X „ _ 

2(e + 7)^1 + 77 (- + - ^^^-^ Si = 



I ,d 1, 



j(i + 1) 



M dr r' 
h 



dr r r 



(59) 



(60) 



(61) 
(62) 



Equation for the big component can be written as Schrodinger equation for a scalar particle 

Si = . (63) 



+ ^ ^ +2M(e+ ") 
dr^ r dr r 



Corresponding non-relativistic 3-dimensional wave function for states with parity P = (— 1)-'^^ is 



.ly+i 



^ = e" 



-lEt/h 



+($i + iSi) S_i 


-(^>i + iSi) S+1 



(64) 
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Now let us consider radial equations for states with opposite parity P 



(4- + -)E2 - 2-Ei = M$o , + -)Ei + + -)Hi = Af$i , 

ar r r r ar r 

+i{e + -)E2 - 2i-Hi = M^2 , + -)^i + -^'o = MEi , 

-i{e + -)^2 - = ME2 , -i{4- + -)^i - = MHi . 

r ar ar r r 



(65) 



Among four dynamical function ^1,^2, Ei,E2 separation of big and small constituents is per- 
formed as follows 



$1 = Si + Ml , ^2 = B2 + M2, iEi = Bi-Mi, iE2 = B2 - M2 ; (66) 
the non-relativistic 3-dimensional wave function for states with parity P = (—1)-' is defined 



(67) 



($1 +i^i) D^i 
P={-iy, xI/ = e-^*M (^^^ + iE,)Do 

{^i+iEi)D+i 

Excluding from (j65p non-dynamical variables <l?o,-ffi, we obtain the system (the rest energy is 
taken away as well: e e -|- M) 

•/ ,^ 1 , d 1. ^ , d 1,^ Z^^n 

i{e + M + -) Ei + —{— + -)[{— + - $1 + -$2 = M$i , 
r M dr r dr r r 

21^ , , d 1,^ I/,, 

t{e + M + -)E2 - xr :r + + -^2 = m^2 , 

r Mr dr r r 

r ,d 2,^ 2i/ ^ , 

-z e + M + - ^>i + — - — + -)E2 El = MEi , 

r Mr dr r r 



C(^^ Id. ,d 2,^ 2u ^ , 

-lie + M + - $2 - T7:r - :r + ^1 = 

r M ar dr r r 



Taking into account (j66|l . transform (j68l) into 



(6 + M + -)(Si - Ml) + + -)2(Bi + Ml) + 

r M ar r 

+17(t- + + M2) = M(i?i + Ml), 

M ar r r 

(e + M + -){B2 - M2) -^{^ + -){Bi + Ml) - 
r Mr ar r 

-^(^2 + M2) = M(52 + M2), 

(e + M + -)(Bi + Ml) - — (^ + -)(S2 - M2) - 
r mr ar r 

,2 



2z/" 

'Mr2 



(Bi - Ml) = M{Bi - Ml] 



(, + M + -)(B, + M,) + --(- + 

+ -)(B2 - M2) + ^~{Bi - Ml) = M(B2 - Ah). 

r M ar r 



(68) 



15 



From whence, we get 

+77(t- + -)-iB2 + M2) = +2M Ml, 
M ar r r 

{e + ^)[B,+M,)-^^{± + l){B,-M,)- 

{Bi - Ml) = -2MMi, 



r Mr ar r 

-M;3(^2 + M2) = +2MM2, 

ie+-)iB. + Ah) + --(- + 

+-){B2 - M2) + ^~{Bi - Ml) = -2M M2. 
r M dr r 

Now summing and subtracting equation within the first couple, and doing the same within second 
couple, we arrive at 

(6 + -)iBi - Ml) + + -)\Bi + Ml) + ^(-f + -)kB2 + M2) + 

r M dr r M dr r r 

+(e + ^)iBi + Ml) - -^(-^ + -)(i?2 - M2) - ^{Bi - Ml) = , 
r Mr ar r Mr"^ 



1 , d 



1, 



, d 



1.1 



(e + -){Bi - Ml) + + -)\Bi + Ml) + + -)-(^2 + M2) - 

r M dr r M dr r r 

-(6 + ^)(i?i + Ml) + -^(-^ + ^)(B2 - M2) + 

r Mr ar r 

+^iBi - Ml) = +4M Ml , 



{e + ^){B2 - ^2) - |^(|: + ^)(i?i + Ml) - |^(i?2 + M2) + 

a,,„ , 1 d , d 2,,„ , 2z^dl,„ 

+(e + -){B2 + M2) + — — {— + -){B2 - M2) + — — -(Si - Ml) = , 



M dr dr 



M dr r 



(6 + -){B2 - M2) - + -){Bi + Ml) 

r Mr dr r 



2u' 



M.2(^2 + M2 

/ ct. , ^ 1 d , d 2.. 

{e+-){B2 + M2)---{- + -){B2 



M dr dr 



M2)- 
+4M M2 . 



Now, using the same method as in 



21 



22l | (consider Bi , B2 as big, and Mi , M2 as small) , we 



arrive at two equations for big components, and two equations defining small components through 



big ones: 



1 , d 1,0 u d 1,1 u d 2, 2i/2 



2iy.d 1,^ 2z^2 1 d d 2, 2vdl^ 

]^^(^ + 7)^1 - J^B, - --(- + -)B2 - J^^-Bi = +4M M, , 



a 2u d 1 



2i/2 1 d , d 2 2u d 



a. 



I , d 1, 



u ,d lA 



V , d 



2v' 



Si = 0. 



Two last equation provides us with non-relativistic radial equations - they can be written 

,2n 



d^ 2d a, 2v 



d^ 2 d , a, 2v'^ 
— + Z— + + -) 2" 



dr"^ r dr 

It is convenient to presents eqs. (I69p in a matrix form 



B2 = 2B2 + AuBi, 
Bi = 2u B2 . 



The right-hand part can be brought to a diagonal form 



Bi 









Bi 


B2 




2v 1 




B2 



fi 




a c 




Bi 




fi 




Ai 




fi 


















f2 




d b 




B2 




f2 




A2 




f2 



The problem is to solve two systems 



a c 
d b 



from whence it follows 



u 
2v 1 



Ai 
A2 



a c 
d b 



Ai a — 2i/ c = 

-u a + {Xi-l) c = , 



Ai 

c = — a ; 
2iy 



\2d-2vb = Q 

-i^ d + (A2 - 1) 6 = 



1- Vl + 4j0- + l) 



-3 , 



b=^d 
2u 



The transformation matrix we need is given by 



fl 




a Ai a/2z^ 




Bi 


f2 




d A2 d/2iy 




B2 



Thus, the system (|70l) is led to the diagonal form 

2d a, 2i/2 2Ai" 



/i=0, 

/2 = 0. 



Note simple relations 



2u 



, 2Ai _ j(i + 1) + 2(j + 1) _ (j + l)(j + 2) 



2z^ 

^2 ^ 2 



2A2 j(i + 1) - 2j (j - i)j 



Thus, we have two problem of one the same type (below to /i and /2 correspond v = j 
V = j + 1 respectively) 



d^ 2d a, + 1) 

— + -- + 2M(e + -)- ^ ^ > 



Changing the variable x = 2\/ —2eM r 



d2 ^2d__l 



aM 



dx^ xdx 4 x^/-2eM 



f = o. 



fix) = 



and introducing the substitution f{x) = x"'e F{x), eq. (j75]) is brought to 



x—r + 2a + 2 - 26x — + 



a(a + 1) - 1/(1/ + 1) 



25 - 2a6 



aM 



V-2eM 



+ 



When 6 = +l/2, a = eq. (f76|) is simplified 



x^^ + {2u + 2-x) 



1 + u 



aM 



F = 



dx"^ ' ' ~ dx y ' ■' yJ-2eM 
what is the confiuent hypergeometric equation for F{A, C; x) with parameters given by 

aM 



A = 1 + u 



V-2eM ' 

The quantization condition is A = — n, which gives 

aM a^M 



C = 2v + 2. 



me 



\/-2eM 



-n 



2{l + v + n) 



2n^il + u + n) 



remembering that to linearly independent solutions with parity P = (—1)-' correspond v 
and V = j + 1. 
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VII. CONCLUSIONS 

Quantum-mechanical system - spin 1 particle in external Coulomb field is studied on the base 
of the matrix Duffin - Kemmer - Petiau formalism with the use of the tetrad technique. With 
the help of parity operator, the radial 10-equation system is divided into two subsystem of 4 and 
6 equations that correspond to parity P = (—1)^+^ and P = {—ly respectively. The system of 
4 equation is reduced to a second order differential equation which coincides with that arising in 
the case of a scalar particle in Coulomb potential. It is shown that the 6-equation system reduces 
to two different differential equations for a "main" function. One main equation reduces to to 
a confluent Heun equation and provides us with energy spectrum. Another main equation is a 
more complex one, and any solutions for it are not constructed. In radial equations, transition to 
non-relativistic case is performed. In this limit, three types of linearly independent solutions have 
been constructed in terms of hypergeometric functions. 
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